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We present cross section predictions for squark and gluino production at the LHC, in association 
with up to two additional hard jets. These cross sections can be very large in comparison to the 
inclusive Born rates. Because hadron collider experiments utilize hard jets in the reconstruction 
of cascade decays or as a way to separate squark and gluino production, the understanding of 
these processes is crucial. We show to what degree hard jet radiation can be described by shower 
algorithms and point out how tuning these showers, for example to top quark pair production, could 
help reduce theoretical uncertainties for new physics searches at the LHC. 



In the Standard Model (SM), the mechanism of the 
observed electroweak symmetry breaking is widely be- 
lieved to involve a Higgs boson. This fundamental scalar 
poses a theoretical problem: the stability of its mass 
after the inclusion of radiative corrections. A possible 
solution is TeV-scale supersymmetry. The minimal su- 
persymmetric extension of the Standard Model (MSSM) 
simultaneously solves several problems in high energy 
physics and cosmology: gauge coupling unification; ra- 
diative electroweak symmetry breaking [ij- and a stable 
weakly- interacting dark matter candidate [^l . 

With the Tevatron in operation and the LHC only a 
few years distant, the TeV scale is rapidly coming within 
reach. Squark and gluino production cross sections ap- 
proach 0{ph) at the Tevatron and 0(nb) at the LHC, 
for masses around the present Tevatron exclusion limits 
of up to 400 GeV. 

MSSM searches and jets: The main difference between 
R-parity-conserving MSSM signals and SM QCD back- 
grounds at a hadron collider comes from the stable light- 
est supersymmetric particle, an end product of all super- 
partner decays, which escapes the detector unobserved. 
Requiring a large amount of missing transverse energy is 
thus the first ingredient to enhance the signal. The QCD- 
strength production channels for squarks and gluinos are 
pp ^ gg, qq* , qq, 0, 13 • To a good approximation, the 
light-flavor MSSM squarks are mass degenerate, while 
the lighter of the two top squarks is often the light- 
est strongly-interacting MSSM particle. Searching for 
squarks and gluinos in an inclusive analysis the signa- 
ture is jets plus ^T, possibly plus leptons. The shortest 
cascade of two- particle decays is g qq and q qx^, 
where the lightest neutralino Xi is stable. Such an inclu- 
sive search is well-suited to find MSSM-type deviations 
from the Standard Model. Because the gluino decay gives 
one more hard jet in the final state, an event's jet multi- 
plicity provides discrimination between the relative rates 
of squarks and gluinos. 

We can also make use of longer decay chains, e.g. of 
the classic type g ^ qq ~* i<^qq ^ tiqq — > Xi^^lQ with 



five unknown masses. These masses can be extracted 
from kinematic distributions, i.e. thresholds and edges of 
different momentum combinations [HI . Alternative meth- 
ods have been developed to improve the mass reconstruc- 
tion [gJ and the associated statistical and systematic er- 
rors. These measurements can in turn be used to deter- 
mine parameters of the TeV-scale MSSM Lagrangian [3] ■ 

To make optimal use of the achievable statistical preci- 
sion, as well as to quantify the errors on extracted model 
parameters, it is crucial to properly understand the sys- 
tematic errors in the cascade reconstruction. Obviously, 
effects such as detector resolution and jet energy scale 
calibration will have a large impact [8]. In this letter, 
we focus on another source of (combinatorial) error: the 
presence of additional, observable hard jets due to SM 
QCD radiation, with transverse momenta comparable to 
the typical cuts planned for squark and gluino studies, 
about pT > 50 - 100 GeV [i, [l3|. Extra jets, in par- 
ticular from the initial state, will introduce noise when 
reconstructing the cascade kinematics from the observed 
jet kinematics and when attempting to separate squark- 
and gluino-enriched samples. This question has a coun- 
terpart in the hadronic top analyses at the Tevatron. We 
know from data that in these analyses additional jet ra- 
diation fakes W decay jets in a non-trivial fraction of the 
events il^. With generally larger energy scales involved 
in SUSY processes, one would expect the initial state to 





CTtot [pb] 


gg ULg ULul ULUL TT 


PT,j > 100 GeV 


a2] 


4.83 5.65 0.286 0.502 1.30 
2.89 2.74 0.136 0.145 0.73 
1.09 0.85 0.049 0.039 0.26 


PTj > 50 GeV 


(Tlj 

a2] 


4.83 5.65 0.286 0.502 1.30 
5.90 5.37 0.283 0.285 1.50 
4.17 3.18 0.179 0.117 1.21 



Table I: Cross sections for the production of a toy-model 
600 GeV top quark, squarks and gluinos at the LHC at the 
benchmark point SPSla. We show fixed-order matrix element 
results with 0,1,2 additional hard jets, and with two different 
values and \yj \ < 5 and Rjj > 0.4. 
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be more active, hence potentially more dangerous. 

Hard jets from matrix elements: To study the produc- 
tion of hard jets with squarks and gluinos at the LHC, 
we first present the results of a fixed-order approach. We 
use the new supersymmetric version [13, of the event 
generator MadEvent [3] to calculate tree-level rates for 
PP ~^ ggA't AiAg^ including the emission of one and 
two additional jets. To avoid regions sensitive to the soft 
and coUinear singularities of initial- and final-state radi- 
ation we limit ourselves to pt,j > 50 GeV (and in some 
cases pT.j > 100 GeV), for which we expect fixed-order 
perturbation theory to be reliable at the LHC. 

In Table U we show inclusive production cross sec- 
tion estimates for squarks and gluinos plus zero to two 
hard jets, for the parameter point SPSla, where rrig — 
608 GeV and m^^ = 567 GeV [l5|]. The factorization 
scale is set to the average final state mass, as is the renor- 
malization scale for the heavy pair. The renormalization 
scale for additional jet radiation is prj, the standard 
choice in shower Monte Carlos We see that for jets 
with pt > 100 GeV the perturbative expansion is stable, 
but the relative suppression is closer to 1/2 • • • 1/3 than 
to as/iT. If we allow semi-hard jets down to 50 GeV the 
fixed-order perturbative expansion approaches its limit. 
We emphasize that it is not a problem for the one-jet rate 
to be slightly larger than the all-inclusive rate at leading 
order, as long as it is in the range of the NLO inclusive 
cross section jj] . To check that this behavior is quantita- 
tively universal for heavy QCD production we also show 
the rates for a heavy toy-quark T (effectively a 600 GeV 
top), and we indeed see the same pattern. 

The case of gluino pair production exemplifies that 
gluon radiation (mostly from the initial state) is indeed 
the dominant source of extra jets: allowing just initial- 
and final-state gluons with pTj > 50 GeV, the gg case 
in Table U becomes CTo,i,2j = (4.63, 3.90, 2.03) pb. In 
addition, jet radiation probes new initial states through 
crossing of quarks and gluons from the final state to 
the initial state: e.g. purely quark-initiated processes 
like qq — > qq receive large corrections from qg scatter- 
ing. These crossed channels include intermediate states 
of the kind qg — > qg — > qqq. Usually these intermediate 
states are subtracted in the narrow-width approximation 
to avoid double counting between NLO qq and gq pro- 
duction Here, we explicitly remove all on-shcU in- 
termediate gluinos from gu^ production. For the ulu)^ 
channel we simply decrease the gluino mass to 558 GeV 
to avoid intermediate on-shell gluinos. 

Jets from parton showers: To compare to experimen- 
tal observables, multiple soft emission and hadronization 
effects can be important - aspects which go beyond the 
scope of a fixed-order calculation. For bremsstrahlung 
emission, logarithms log^^('9soft /'Shard) appear to 
all orders in perturbation theory, while the correc- 



tions associated with hadronization are inherently non- 
perturbative, — > 1. Exploiting coUinear factorization 
in QCD and assuming universality of hadronization, both 
types of corrections may be included, at least approxi- 
mately. This is the basis of multi-purpose event genera- 
tors like Herwig, Pythia, or Sherpa. Starting from a 
given set of final state partons, a sequence of initial- and 
final-state QCD branchings are generated, resumming 
the leading logarithms mentioned above. The emissions 
are ordered, e.g. in the parton virtuality Q or in pT,j^ and 
the description is matched to hadronization models at a 
fixed low scale, '-^ 1 GeV. 

Especially for hard radiation, large differences may ex- 
ist between different shower algorithms. To quantify, we 
use Pythia 17[ with two qualitatively different show- 
ers: one Q^-ordered and the other pT-ordered 
Note that, due to the large final state squark and gluino 
masses, we explore mainly the properties of the initial- 
state showers. The crucial parameter here is the starting 
scale of the shower, which sets an upper limit to the phase 
space over which jets can be radiated. For initial-state 
radiation, this starting scale is nominally identical to the 
factorization scale, where the parton densities are convo- 
luted with the matrix elements. In Pythia, this scale is 
normally the transverse mass fip = \/Pt '^ith 
the average mass of the final state SUSY particles, and 
Pt their relative transverse momentum. 

For a pr-ordered shower, this ixp can be used directly 
as the maximum ptj - We refer to this choice as the 
PT-ordered 'wimpy shower'. We also investigate the con- 
sequences of allowing the parton shower to populate the 
full phase space, with the maximum p^.j = -s/s/^i re- 
gardless of /XF. This choice we refer to as the pr-ordered 
'power shower'. Though strictly-speaking in conflict with 
the factorization assumption, this choice has interesting 
phenomenological consequences, as we shall see. 

The case of a Q^-ordered shower is not so simple. The 
starting scale here is Q^ax = {Cf^p,s) , where C > 1 
parameterizes the translation from p|, to Q^. We refer 
to C = 1 as the Q^-ordered wimpy shower, C = 4 as 
Tune A (designed to match Tevatron data) , and C 
oo as power shower — with the same caveat concerning 
factorization as for the p^-ordered version. 

Numerical Comparison: We now turn to a comparison 
of the three processes tt-|-jets, g^-l-jets, and UiUi-l-jets, 
at the LHC. Here, tt represents a process with a small 
ratio fiF/\/s, while the SUSY processes involve much 
larger masses and hence larger factorization scales. The 
difference between gg and ulul is the sea- v. valence- 
dominated initial state. 

Jets in the parton shower final states were clustered us- 
ing a cone algorithm with Ai? — 0.4, i.e. similar to the 
ME cut Rjj > 0.4. By looking at ptj/pt.me in gener- 
ated dijet samples, we checked explicitly that out-of-cone 
corrections are too small to affect our study significantly. 
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Figure 1: prj spectra for tt,gg and ulul production in association with 1 and 2 hard jets. NLO K- factors [J] are applied 
to the Pythia resuhs to avoid a normalization mismatch in the one-jet case. The mass spectrum is given by SPSla, except 
for ulul where we reduced the gluino mass, as in Tab.|ll At high pr (of order the factorization scale and above), the matrix 
elements are, by definition, the most reliable while the parton showers can be seen to be associated with large uncertainties, 
whereas the opposite is the case at low pr- The shaded region in the middle row is the theoretical uncertainty on the matrix 
element prediction. See text for details. 



In Fig. [T] we show results for the pT.jspectra, using 



CTEQ5L parton distribution functions |2l|. We use the 
heavy particle mass as the central value for the factoriza- 
tfon scale. For renormalization scale, we apply one factor 
of as using — m3p for each final-state particle. The 
shaded band in the middle row highlights the theoreti- 
cal uncertainty by plotting the extremal values of four 
factorization and renormalization scale choices: varying 
pLR down to the minimum pT of the jets, and varying 
down to the minimum pT of the jets and up to \/i/2. 

Consider first tt-|-ljet (upper left): as a general fea- 
ture, the power showers exhibit a harder high-p^ spec- 
trum than the matrix element. In contrast, note the 
cataclysmic drop of the wimpy showers above pT,j 
"fTit ^ ^jLf, with the Q^-ordered wimpy showers every- 
where softer than the matrix element. Finally, Tune A 
is indistinguishable from the matrix element for pTj < 
250 GeV, this region being similar to the phase space ac- 
cessible at the Tevatron. All the calculations agree fairly 
well down to pxj ~ 50 GeV, a value sufficiently below the 
factorization scale for the parton showers to be reliable, 
but still large enough for stable fixed-order predictions. 
The crossover at 150 GeV between the two pT-ordered 
showers illustrates the effect of changing from pt/2 
in the wimpy shower to 3pT in the power shower. The 

-ordered showers all use = PT,j- 



For the two-jet distributions, note that the coUinear 
approximation is only rigorously correct in the limit that 
each successive jet is much softer than the preceding 
one, hence the two-jet shower predictions are associated 
with even larger uncertainties. Nevertheless, the power 
showers deviate from the matrix element by less than 
a factor two over most of the pTj range. Although re- 
maining differences could undoubtedly be cured by slight 
shower parameter modifications, we note that obtaining 
simultaneous agreement of the one- and two-jet rates 
is likely to require more sophisticated matching tech- 
niques [Ullilii,!!?!. 

The parton shower predictions for the SUSY processes 
exhibit significantly less variation, as shown in the lower 
rows of Fig. [T] Owing to much larger factorization scales 
the presence or absence of a cutoff in the parton shower 
evolution at /if does not lead to very large differences for 
the pT,j regions we consider. In this way, the kinematic 
regime of squark and gluino production at the LHC is 
more similar to e.g. tt production at the Tevatron than 
at the LHC. Considering first gg+1 jet, we observe that 
the matrix element rate begins to diverge around px^j ^ 
100 GeV, due to large logarithms, log^m|/py j. For Z, 
light Higgs, or even tt production this breakdown scale 
is much smaller, of order the minimum pxj observable 
at the LHC [2^ . The case of heavy particle production. 
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Figure 2: The A/Jjj distribution for ggjj production as pre- 
dicted by the hard matrix element and the parton showers. 

such as gluinos and squarks, is different: for jets softer 
than pt ^ 100 GeV we would be well-advised to include 
resummation effects. 

At large values of pxj , a pattern similar to the tops 
arises: the one-jet radiation is very well described by 
Tune A. For two-jet radiation, the Q^-ordered showers 
generally fall below the matrix element, while the px- 
ordered ones overshoot. We also studied mlu^ and ULg 
production and found that they exhibit essentially the 
same behavior as gluino-pair production. 

For the last process in Fig. [TJ Ulul production, the 
matrix element divergence at low px appears to be much 
milder than for the gg case. We interpret this as a 
consequence of the less radiating, valence-dominated ini- 
tial state. For the high-pT tail, the power showers are 
again in fairly good agreement with the matrix element, 
though with a much smaller difference between px- and 
Q^-ordered showers than above. 

Lastly, we compare the A.Rjj distributions for gluino 
pairs plus two jets in Fig. [21 The results are not dras- 
tically different. The difference between the two differ- 
ent shower models at low AR is interesting, however. 
It should probably not be interpreted as the onset of a 
coUinear singularity. In that case, one would expect the 
showers to agree with each other, but not with the ma- 
trix elements. Moreover, from the results presented in 
Fig. [1] it is clear that with a jet cut of 100 GeV we are 
nowhere near the coUinear region. In fact, the cutoff it- 
self may furnish part of the reason; the region where both 
jets are close to the cutoff is, by definition, not strongly 
ordered, and would hence be expected to be problematic 
for shower descriptions. We plan to return to this in a 
future study. 

Summary: Using SUSY-MadEvent we show that 
matrix-element-based QCD calculations predict a large 
number of hard jets associated with squark and gluino 
production at the LHC. This effect should be taken into 
account in studies of the separation of squark and gluino 
event samples, and for cascade decay reconstruction. We 
compared in detail the matrix element approach and the 
results from p^-ordered and Q'^-ordered showers, imple- 
mented in Pythia 6.3. 

For the radiation of one extra jet in SUSY heavy 
colored pair production processes, conventional parton 



showers with a phase space cutoff at the factorization 
scale give a reasonable approximation up to pT,j ~ IJ'f/'^, 
above which they rapidly break down. This is similar to 
what has been observed before for Drell-Yan 
(and Higgs) production. A significant improvement can 
be obtained by removing the explicit phase space cut, 
allowing the shower to populate all of phase space. How- 
ever, this tends to yield somewhat harder radiation spec- 
tra than produced by the matrix elements, again akin 
to what has been shown for hadroproduction of colorless 
resonances [291] . 

For fairly soft jets, we see that in the production of 
high-mass gluinos the breakdown of fixed-order perturba- 
tion theory caused by logarithmic corrections can occur 
already at jet transverse momenta of as high as 100 GeV. 
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